Atomic Josephson vortex 
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We show that Josephson vortices in a quasi-lD atomic Bose Josephson junction can be controUably 
manipulated by imposing a difference of chemical potentials on the atomic BEC waveguides forming 
the junction. This effect, which has its origin in the Berry phase structure of a vortex, turns out to 
be very robust in the whole range of the parameters where such vortices can exist. We also propose 
that a Josephson vortex can be created by the phase imprinting technique and can be identified by 
a specific tangential feature in the interference picture produced by expanding clouds released from 
the waveguides. 
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I. INTRODUCTION 

From the very beginning of the experimental achieve- 
ment of Bose-Einstein condensation in trapped ultracold 
gases, various solitonic and topological collective struc- 
tures, such as dark and bright j2| solitons, vortices 
(3, skyrmions 0|, etc., became objects of great theoret- 
ical and experimental interest. These configurations are 
direct manifestations of macroscopic coherence of BEC 
and many of them were repeatedly created and detected. 

In our previous work Q it has been shown that yet an- 
other stable object, namely atomic Bose Josephson vor- 
tex (JV), can exist in a long Bose Josephson junction 
(BJJ). Such junction can be formed between two parallel 
quasi- ID Bose Einstein condensates coupled by tunnel- 
ing. The circulating atomic supercurrents — counter- 
propagating in each waveguide and closing the loop be- 
tween them — represent a neutral analog of the Joseph- 
son vortex studied in superconducting Josephson junc- 
tions in great detail While bearing similarities with 
the charged case, the absence of the Meissner effect makes 
the local formulation of the BJJ solution possible only in 
a quasi-lD geometry. The most interesting feature of the 
Bose JV is its interconversion into a dark soliton (DS) ■ 
Accordingly, the description of the Bose JV necessarily 
involves both the phase and the amplitude of the wave- 
function. This instability, which can be controlled by 
tuning the Josephson coupling, may potentially be uti- 
lized in several ways 5]. Here we show that the Berry 
phase term is responsible for a force produced on the JV 
by a time-dependent difference of the chemical potentials 
applied between the waveguides. Once the JV speed 
as a whole reaches a certain critical value, the instability 
destroying the supercurrent circulation develops so that 
the JV transforms into a moving DS (grey soliton) which 
is essentially insensitive to Sfi. We also discuss how the 
phase imprinting method and the interference after some 
free expansion can be employed for creating and detect- 
ing the JV. As it turns out, the formation of the relative 
phase between the waveguides is quite insensitive to the 
characteristics of the imprinting beams once certain topo- 
logical requirements are satisfied. The JV introduces a 



specific tangential feature into the interference pattern, 
which should allow an unambigious identification of the 
JV as well as of the interconversion effect [5| . 

Experimental tools with great potential for creation, 
manipulation and detection of the JV have already be- 
come available and continue to develop. Long BJJ con- 
sisting of two parallel coupled atomic waveguides was ex- 
perimentally realized with two different and equally ver- 
satile techniques capable of producing double- well poten- 
tial of necessary geometry. First technique uses counter- 
propagating laser beams to form an optical dipole trap 
!(| and second one uses a microfabricated atom chip with 
purely magnetic potential In both cases relative co- 
herence of BECs has been revealed in the interference 
fringes. 



II. ATOMIC QUASI-ID JOSEPHSON VORTEX 

Below we will introduce a phenomenological descrip- 
tion of the atomic JV which captures three main effects: 
i) The JV instability with respect to its interconversion 
into the DS as the Josephson coupling strength exceeds 
certain threshold |^; ii) The critical speed effect; iii) The 
effect of force on the JV produced by uniform We 
will also derive the corresponding phenomenological co- 
efficients in terms of the microscopic model Q close to 
the threshold and will consider the opposite limit — very 
small couplings, when the JV can be well described by 
the Sine-Gordon (SG) equation. 



A. Phenomenological free energy 

Phenomenological description helps elucidating generic 
features which do not depend qualitatively on a particu- 
lar microscopic model. The main characteristic of the JV 
is a presence of a localized supercurrent circulation de- 
generate with respect to its orientation. This circulation 
can be decsribed by persistent currents Ji^2{x,t) along 
the first and the second waveguides, respectively, as well 
as by tunneling Josephson currents between the waveg- 
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uides. Here a; is a spatial coordinate along two parallel 
waveguides and t stands for time. If the JV is stationary, 
the currents Ji.2{x,t) are localized within some typical 
distance and flow in opposite directions, and, therefore, 
the difference J{t) — J dx{Ji{x,t) — J2{x,t)) is a good 
global measure of the circulation. Obviously, the quan- 
tity P{t) = J dx{Ji{x,t) + J2{x,t)) is linear momentum 
of the JV. It is zero, if the JV is stationary, and it must 
acquire a finite value, if it is moving as a whole. The 
total momentum is assigned to center of mass of the JV 
positioned somewhere at XQ(t) along the junction. This 
generic picture is well known from the phenomenological 
description of the JV in superconductors (8|, where the 
SG model is sufficient to find the relation between P{t) 
and the center of mass velocity V = Xq through the JV 
effective mass M. 

It is important that in the SG model the circulation is 
essentially a topologically conserved quantity (provided 
the contribution to the circulation from the tunneling 
currents are negligible). Thus, this model does not con- 
tain the intercoversion effect p| taking place in a quasi- 
ID geometry. Incorporating density fluctuations shows 
that the center of mass position Xq can be associated 
with a dip in the condensate densities in each waveguide. 
This dip becomes more pronounced as the Josephson cou- 
pling 7 strengthens, and, finally, it becomes zero of the 
densities at some critical value 7 = 7c, which signals a 
moment of the transformation of the JV into the DS. At 
this point and for higher 7 no more current circulation J 
exists. Obviously, this effect can be interpreted in terms 
of spontaneous breaking of time reversal symmetry for 
7 < 7c with the order parameter J ^ 0. Since J is just 
one degree of freedom the corresponding "phase tran- 
sition" occurs in zero spatial dimensions and therefore 
cannot be considered as a true phase transition. Yet, the 
corresponding free energy of the system containing just 
one DS or JV can be represented in a sense of the Landau 
expansion of the effective action close to the "critical" 
point 7 ~ 7c as 



^eff = XqP + QJ — Tieff, 
with the effective Hamiltonian being 



(1) 



p2 



2M+ 2M_ 



-+a{-f--fc)J^+ClJ'^-C2SllPJ+C3J^P^. 

(2) 

where Q{t) stands for a variable canonically conjugate 
J; M+, Af_, a > 0, ci > 0, C2, C3 > are phenomeno- 
logical coefficients dependent on a particular form of a 
microscopic description. These coefficients will be deter- 
mined later within the variational approach applied to 
the model 0. 

The Lagrangian I|1I2|I describes the effect i), that is, 
spontaneous formation of the circulation J 7^ for 7 < 7c 
I3 (ci > insures stability beyond linear approximation) 
as well as coupling between the center of mass motion and 
the internal circulation. 

The term ~ C3 accounts for the effect ii): a.s P ^ V 
exceeds some critical value, the effective coefficient 0(7 — 



7c) + caP^ becomes positive, which restores the time- 
reversal symmetry so that the circulation J — 0. 

The term ~ C2 is responsible for the effect iii), that is, 
for the force induced by Sfi. This term is symmetric with 
respect to exchanging of the waveguides and conforms 
with the time-reversal symmetry. Indeed, swapping the 
waveguides positions changes 5/i — > — 5/i and J — > — J. 
The time-reversal leaves the product JP invariant. The 
term SfiPJ is responsible for the force ~ 6fiJ on the JV. 
Its origin can be understood as follows: When 0, 
a number of atoms ~ Sfj, starts tunneling between the 
waveguides. Accordingly, the first and the second waveg- 
uides acquire the total linear momenta ~ (5/i / dxJi and 
^ —dfi J dxJ2, respectively. Hence, the total momentum 
P attains a nonzero value ~ SfiJ as long as the circula- 
tion J is finite. The force, then, is P ^ 6fiJ. It can also 
be viewed as a force due to the bias introduced by the 
tunneling currents driven by (compare with the bias 
term in the SG equation 8J). 

In the presence of dissipation one should introduce the 
dissipative function in terms of the velocities Xq and Q 
as 



2a+ 



2<tJ 



(3) 



with some kinetic coefficients a± > Q. Then the stan- 
dard variational procedure with respect to the conjugate 
variables yields the equations of motion 

V-(^-^ + 2c3J^^P-C25^iJ^0, (4) 

P+ — =0, (5) 

Q - 2 (a(7 - 7c) + J - 4ci J3 - C26fiP = 0, (6) 

j+^ + A^O. (7) 

This system describes a quite complex non-linear dynam- 
ics of the circulation J coupled to the center of mass mo- 
tion. Besides small oscillations of J around equilibrium, 
strongly non-linear evolution of J corresponding to the 
phase-slip, during which J changes sign, can occur as 7 
approaches 7c from below. The term ~ C3 in eq.|0 is 
responsible for strong coupling of the phase-slip to the 
motion of the JV as a whole. Furthermore, the terms 
~ 6fj, introduces the parametrical driving. Thus, all three 
effects mentioned above cannot in general be well sepa- 
rated from each other, and general solution of the system 
can only be obtained numerically. However, considering 
certain limiting cases is helpfuU for making these effects 
independent from each other. So, in what follows we will 
consider such limiting situations. In reality, of course, 
the above system should be analyzed numerically. 

Let us, first, consider the case S/j, = and P — 0, 
that is, stationary JV. Then, eas. H6l7|l are essentially 
those discussed in ref.j5j] which describe the dynamics 
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of the JV-DS interconversion which involves the phase- 
slip. Here we will not discuss this effect (that is, i)) any 
further. 

The effect ii) makes sense to consider only in the 
case of small dissipation with respect the center of mass 
motion. So, we set \/<j+ = and, keeping 5^ = 
0, find that the generalized momentum is conserved 
P = const. Then, eq.lgj gives P = M+V in the 
limit of small J. Further substitution into eq.© yields 
-2 (a(7 - 7c) + CaMlV'^) J - 4ci = in equilibrium. 
A non-trivial solution exists only when the coefhcient 
a(7 — 7c) -I- c^M^V^ in front the linear term in eq.® 
is negative. Accordingly, the critical velocity becomes 



(C3M2)^ 



(8) 



above which no quasi-static solution for the circulation 
J 7^ exists. 

Now, let's consider the effect iii) where (5/i 7^ and 
small. Then, we assume that velocities are much smaller 
then the critical one (|SJ) and the JV circulation pa- 
ram eter J is in equilibrium close to its value Jq = 
i\/ 0^(70 — 7)/2ci determined when = 0. In this situa- 
tion, one should ignore the term ~ ci in eq. because it 
describes just a small correction to Jq. Subsequently, J 
can be replaced by Jq in cq. 10} and term ~ C3 should be 
dropped because it describes a small correction ~ 7c — 7 
to the finite term 1/M+. Then, excluding P, one finds 
that the center of mass velocity obeys the equation 



V 



AU<T, 



m 

M+ ' 



fit) - C2M+J6fi. 



(9) 



As discussed above, the term f{t) describes the force in- 
duced by time dependence of the difference of the chem- 
ical potentials J/Lt. The relation of this force to the Berry 
phase term in the full "microscopic" action will be con- 
sidered below. 



B. Variational approach 

It is worth noting that all the phenomenological coef- 
ficients can be derived from a "microscopic" Lagrangian 
and a dissipative function. Here we will use a simpli- 
fied approach which ignores dissipation and will consider 
model [3 as the "microscopic" Lagrangian. In terms of 
the fields "01.2 describing each waveguide the Lagrangian, 



is given by the Berry term 

and by the Hamiltonian 

Ti. = Hi + H2 + Ti-u, 



(10) 



(11) 



(12) 



consisting of the each waveguide terms 

■Hk^ [ dx[^\V^,\' + f - (13) 
J 2m 2 

with k ~ 1,2, and of the contribution which describes 
the Josephson tunneling between the waveguides 



H12 = dxj{iplip2 + tp2'^l), 



(14) 



where /ij, are waveguide's chemical potentials and the 
integration J dx... is performed along the waveguides. 

The equations of motion following from eas. H10ll4|l in 
the units h — 1, m = 1, with the imit of length given by 
the healing length Ic = 1/v^ for /i = (/ii + /i2)/2 and 
the unit of time determined by to = l//i, are 



iiPi = — —^i-{i + vU'i + \'i/jif^i-i^ip2; (15) 

iip2 = — - (1 - '7)V'2 + |V'2|^V'2 - J^V'i- (16) 

Here the wave functions were transformed as ipk — > 
^/noi!k, where no — [ijg is the average ID density of a 
single uncoupled (7 = 0) waveguide; the quantity v repre- 
sents the dimensionless Josephson coupling v = 7//i, and 
r\ stands for the difference of chemical potentials rescaled 
by 2/^. These equations admit exact JV stationary solu- 
tion discussed in detail in ref.*^. For small velocity of 
the JV, it is natural to use the variational ansatz which 
coincides with the stationary solution at y = 0. Thus, 
we choose 



^-1,2 = 0^ - ^O(O)) + 



i\fsQ 



1,2 



ch(s {x-Xo{t))) 

(17) 

with the parameter s giving the JV size. The bulk den- 
sities ni^2 can be obtained in the thermodynamical limit 
(when no JV is present) from eas. H15ll6(l . Since we are 
interested in small deviations only, the corresponding ex- 
plicit expressions are 



?T.1,2 



(1 



1± 



■n 



1 + 2v 



(18) 



where the indexes 1, 2 correspond to ±, respectively. 

The stationary solution [^j can be obtained from the 
ansatz (|17|l by setting Xq = const, y/sQi = —\fsQ2 = 
±\/l — "iv and s = 2-^ (compare with ref.0]). Con- 
sidering complex Qi,2 and real Xq as slow dynamical 
variables, one can substitute the ansatz (|17|l into the La- 
grangian (|10I14() and perform explicit integration over 
X . This procedure generates the effective Lagrangian 
£e in terms of the variables Qi,2, ^0, s and their time- 
derivatives. Obviously, such procedure is in line with 
separation of fast and slow variables, so that only slow 
dynamics should be considered to full extent. Close to 
the interconversion instability (y k. Vc = 7c//^), the slow 
variables are Xq and J . The variable s describes fast 
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adjustment of the JV size. It is important that it is not 
dynamical within the chosen ansatz. Indeed, as can be 
seen, the effective Lagrangian does not contain s. The 
choice of s is dictated by conservation of total number 
of particles during dynamical evolution of the other pa- 
rameters. Calculating the depletion SN of the number of 
particles caused by the presence of the JV, we find 



SN 



2(ni + 712) 



2(|QiP + |Q2n=C 



(19) 



where the constant Cat is determined for the station- 
ary JV by setting all the time derivatives to zero and 
minimizing the effective energy with respect to Qi^2 
and s. Considering small values 77, it is enough to set 
ni + n2 — 2(1 -|- v), which is the equilibrium value. Here 
we will consider values Qi,2 — *■ 0, so that the explicit 
solution of ea. lfTO|) for s becomes 



1 - 3j/ 



2^/T 



IQ1P + IQ2 



(20) 



As discussed in ref.^, the value = = 1/3 is the crit- 
ical point below which the JV forms spontaneously from 
the DS. Thus, the smallness of Qi,2 automatically implies 
a proximity to the critical point. Then, for consistency 
of the effective action expressed in powers of Qi,2, the 
value v should be set to Vc except in the quadratic term 
vanishing at the critical point. 

It is worth discussing, first, the structure of the Berry- 
term part l|ll|) of the full action. As mentioned above, 
the cross term ~ r/PJ leading to the force on the JV 
~ ?7 in the Lagrangian |^ can be viewed as generated 
by the Berry phase effect. Indeed, the Berry part is 
Cb = J dx{—pi(pi — p2^2), where pk and ipk are den- 
sity and the phase, respectively, in the fc-th waveguide. 
In the static solution as well as in the ansatz H17|l 
each phase changes by ±7r, so that, e.g., if at a; = —00 
one finds ipi = ip2 = 0, then, at x = +00 there 
is ^fi — = Svr. If the JV is moving slowly, then 
ifk ~ XoV</Jfe. Thus, a substitution into the Berry part 
gives Cb XqJ dx{piV(pi + p2^'~P2) « 7rXo(pi - P2), 
where spatial variations of the densities are ignored. Flip- 
ping the time derivative and realizing that pi — P2 ^ ^, 
one finds Cb ~ "A^o'ij which is the work done while 
making a displacement by the force ~ 77. 

The relation between the observables P, J can be ob- 
tained as a result of substituting the ansatz (|17|l into 
(inj- This yields Cb = XqP - 2{B+Q+ + B^Q^), 
where Q+ = Re(Qi -I- Q2), Q- ^ Re(Qi - Q2), B+ = 
Im((5i + Q2) and B- — lm{Qi — Q2)- The total momen- 
tum P = ~i J dx{iplT/'ipi+'ip2^'ip2) and the supercurrent 
circulation J = —i J dx{iplV ipi — ^2^^2) are given as 



P = -^{l + vf/^ 



Q- 



vQ- 



2(1 + 2i/) 

vQ+ 



2(1 + 2iy) 



(21) 
(22) 



Thus, the parameters Q± can uniquely be expressed in 
terms of P, J. In particular, ~ P for 77 = 0. The 
Berry phase effect discussed above is reflected in the part 
~ 77 of ea. l|21|) . It is also clear that the variable conju- 
gated to J in eq.(^ is Q ~ B-. Further analysis shows 
that the quantity P-|_ = Im((5i + Q2) enters Cb in a 
combination ^ P+P — 7'bP+ + o(P+) with some > 0. 
Thus, P-i- would generate higher time derivatives with 
respect to Xq, which should be neglected as long as the 
JV motion is slow. Accordingly, it is reasonable to set 
P-(- = in ea.H17|l. so that Qi^2 are chosen in the form 



Q+±Q- ^ .P- 
^^•^ 2 2 



(23) 



with the indexes 1, 2 corresponding to ±, respectively. Fi- 
nally, employing the ansatz (jl7l23(l in the "microscopic" 
Lagrangian H1UI14|I and expressing the variables (|23|l in 
terms of the observables P, Xq, Q, J as described above 
one arrives at the effective Lagrangian in the form 
where the coefficients (in the chosen units) are 



1 

M+ 


V3 
27r2' 


(24) 


1 

mT 


327^2 

27 V^' 


(25) 


-7c) = 


3V3 / 1\ 
4.2 gj' 


(26) 


Cl = 


3V3 
1287r4' 


(27) 


C2 = 


2V^ 
57r2 ' 


(28) 


C3 = 


9\/3 
647r4' 


(29) 



In these expressions, the difference ^ ly ~ i/c has been 
ignored except in the quadratic coefficient (|25|) . which 
determines the instability. 

It is worth noting that far from the instability (that is, 
^ T^c) no simple expansion for the effective Lagrangian 
in terms of powers of Qi,2 can be obtained. However, it 
is clear that in this limit, the variations of the density 
can be ignored as long as the JV is close to its internal 
equilibrium. Thus, the effective action can simply be ob- 
tained by considering fluctuations of the densities being 
small. This will lead to the SG model. For completness, 
we will consider this limit below. 



C. Sine-Gordon approximation 



Setting Tpi 
\/n2 + n'2{x,t) exp(z(/j) 



= \/ni+ n[{x,t) exp{iip), ■02 = 
with ni^2 being uniform back- 



ground densities given by ea. lfl^ and 



n,2 



describing 



small fluctuations, and substituting this into the action 
(|10|l one obtains the effective Lagrangian in terms of the 
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relative phase Lp only as 

SG = J dx[{n2 - ni)tp + tp'^ 

- (1 + J/)(V(^)2 +2z/(l + i^)cos(2(/7)], (30) 



250 



where n'l 2 were eliminated in the long wave limit, with 
the gradients of n[ 2 ignored; the difference n2 — ni 7^ 
has only been retained in the Berry term. 
The variation yields 

(/3 - (1 + z/)V^((9 + 4t/(l + ly) sin(2(^) ^ hi - n2 (31) 

the biased SG equation with its well known solution in 
the static limit (see in *8^). It is important to note that 
this equation exhibits Lorcntz invariance (for hi ~ 712 — 
0). Hence, no the critical velocity effect ii) can be de- 
scribed in the SG approximation. 

Here we will be interested the limit of slow velocity 
V and small hi — h2 in ea. H31|l . Then, one can use the 
static SG solution, with its center of mass velocity V 
being a slow variable in the Lagrangian 1)30(1 . Then, 
integrating explicitly, one finds 



-SG 



J dx[-2TrT]V 



MsgV\ 



where the effective mass is 



Msg = 4:V2v. 



(32) 



(33) 



Varying with respect to Xq and including the dissi- 
pative function, one arrives at the equation similar to 
eq.®: 



V 



(T+MsG 



27rr) 



SG 



(34) 



which describes the effect of force induced by the Berry 
term. Thus, while a particular expression for the force 
depends on the proximity to Vc, the effect persists for all 
values V < Vr. 



D. Numerical simulation of the Berry phase effect 

Sensitivity to the difference of the chemical potentials 
•q can be a useful tool for manipulation of the JV posi- 
tion in the junction. To demonstrate this numerically we 
used T] as an externally controlled variable to displace the 
JV on a distance much greater than its size and, then, 
to return it to its original position. The result of the 
simulations of the full system (|15l I16|l is displayed on 
FIG^ The plot on the left represents density of a sin- 
gle waveguide by the intensity of white color. The dark 
curve is a trajectory of the JV center, where its density is 
minimal. Shown on the right, the time-dependence 77 is 
chosen as ri{t) = O.lsin(wi), with ui = 7r/125. The sim- 
ulations have beenperformed with the dissipative term 
introduced in ref.'5|. It is important to note that the 
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FIG. 1: Motion of the JV along the junction (left) gener- 
ated by change of relative chemical potential t] (right) from 
numerical simulations of the system 11511161 . Left figure: the 
atomic density in a single waveguide is represented by inten- 
sity of white color; the JV path in x-t coordinates corresponds 
to the dark color (the density depletion in the JV center) rel- 
ative to the white color of the background. Here v = 0.1, 
and the dissipative term from ref. 5] is given by the value of 
the kinetic parameter a = 2. The units of x and t are the 
coherence length 1^ and the coherence time to, respectively. 



regime chosen for numerical simulations is neither close 
to the critical point nor to the SG limit. Yet, the time 
dependence Xo(t) = ^(1 — cos((jjt)), which follows from 
eas. (|9l34f) in the limit of small damping and where A is 
some amplitude, is consistent with the simulations with 
the reservation that the numerical value of A is not repro- 
duced correctly by either limiting approximation: while 
eq.® underestimates A by a factor of 2, the SG limit 
(jSljl overestimates it by a factor of 4. 



III. 



CREATION BY PHASE IMPRINTING 



The JV can be formed as a result of the decay of the 
DS, once the Josephson coupling 7 is reduced below a 
critical value 7c (or, in the chosen units, with the crit- 
ical value being = 1/3) 0. An alternative method 
is the phase imprinting. It is already a well established 
experimental tool for creation of the DS 0]. It consists 
of exposing a BEG to a pulse of a far detuned laser beam 
which acts as a temporary external potential U{x,t). Ac- 
cording to the impulse approximation, the duration of the 
pulse 6t must be short compared to the correlation time 
of the condensate to = 1/^ , so that no change of the 
BEG density occurs during the pulse — atoms just ac- 
quire finite speeds without performing any significant dis- 
placements. The phase, on the other hand, accumulates 
according to 5ip{x) = J dtU{x,t) and the wave function 
tp before the pulse is transformed as -0 ^ e~^^^4' after 
the pulse. To create a DS, one needs to expose a half- 
plane of an elongated BEG to a laser pulse with a spatial 
variation reminiscent of the typical DS phase profile (the 
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IV. 



DETECTION BY INTERFERENCE 



FIG. 2: Evolution of the relative phase after the imprinting 
as follows from eas. ll5ll6ll . Small dissipative term discussed 
in ref.vSj with a — 6 has been added; i/ — 0.1. The units of x 
and t are the coherence length Ic and the coherence time to 
respectively. 



TT-step). In order to produce the JV, one needs to apply 
a pulse with spatial profiles Ui^2{x,t) specific for each 
waveguide. These should reflect the structure of the JV 
phases in each waveguide: ipi,2 = at ,e.g., x — —oo and 
ipi = —ip2 = ±7r at x = +0O, with smooth transition in 
between at a typical length comparable to the JV size 
= l/{2y/i'). Accordingly, Ui^2{x,t) = at x = —oo and 
Ui{x,t) = —U2{x,t) at a; = +oo, with the "crossover" 
region being approximately equal to the JV size and the 
time integral J dtUi{x, t) — n at x +oo. 

It is important to note that, once the above "topo- 
logical" requirements are satisfied, the JV forms with 
minimal disturbances regardless of other details of the 
laser beams profiles. The most robust characteristic of 
the evolving solution turns out to be the phase differ- 
ence between the waveguides. In the case of the DS, the 
complete density depletion must form at the DS center. 
Thus, the adjustment is accompanied by a strong pertur- 
bation in the form of the density waves [J . The depletion 
at the JV center is also strong, if i/ is close to the critical 
value 1/3. Accordingly, the densities in each waveguide 
will experience significant perturbations. Yet, the phase 
difference relaxes quite smoothly to the equilibrium pro- 
file. To demonstrate this feature, we ran numerical sim- 
ulations with the initially imprinted profile of the phases 
given by the tanh-type variations of the light intensities 
as described above. The result of the following evolution 
of the phase difference is represented on FIG|21 As can 
be seen, the equilibrium phase profile establishes after 
few relatively small oscillations even though the initial 
extension of the (imprinted) phase was about two times 
longer than the equilibrium JV size. 



Experimental visualization is typically done by absorp- 
tion imaging |^ ^3 1 with its intensity proportional to 
the density n of the expanding cloud. In contrast to bulk 
vortices, which can be detected by observing their cores, 
the JV does not have a core. Yet, the phases exhibit the 
TT-jumps. Thus, the interference of the expanding clouds 
released from the waveguides should demonstrate a cor- 
responding feature. 

In quasi-lD regime a good approximation for the 
waveguides wave functions in transverse directions is the 
Gaussians G{y, z) — exp (— (y^ -I- z^)/2(i^), with d being 
a typical width of each waveguide. Thus, in 3D, the two 
waveguides separated by a distance 2zq can be described 
by the following ansatz: 



*o(i?) 



= f{x)i!i.2{x)G{y,z±zo) 



(35) 
(36) 



where ^'1,2(2;) are the solutions corresponding to 
either the DS {Qi,2 — 0) or to the JV and the sign 
± is different for different waveguides. The envelope 
f{x) = (1 — (4a;^/L)), with L, the axial system size, be- 
ing much larger than the JV, reflects finiteness of the 
BEG clouds. 

As long as the transverse dimension d is much smaller 
than any axial feature, the expansion occurs primar- 
ily in the transverse direction. Thus, the density de- 
creases, practically instantaneously, so that the expan- 
sion is essentially free of interaction. This can be for- 
mulated as the requirement {d/lcY ^ 1. Indeed, the 
mean-field interaction is given by the chemical poten- 
tial /i ^ ?^(^)I where n{t) ^ 1/R^{t) stands for a typi- 
cal density of the expanding cloud scaled by its radius 
R{t) w ^/(P + (J/d)^ (in chosen units). The interaction- 
induced additional phase shift can be estimated as A(j) « 
lidt ~ ndt fa {d/lcY ^ 1; where the initial den- 
sity is taken as n = 1 in the chosen units. Hence, under 
this condition, the density after time t of free expansion 
becomes 



n{R,t) 



1 



d^R' exp 



i{R-R') 
2t 



t\2 



*o(i?') 



(37) 



A numerical factor in is set to 1, because it defines 
only overall intensity (not the structure) of the absorp- 
tion image. 

When two expanding uniform BECs overlap they form 
an interference pattern (IP) of parallel fringes The 
specific signature of a rotational vortex in the IP is a so- 
called edge dislocation. It was predicted in ref.Tl'l and 
then seen experimentally in rcf. 10] . Here we will discuss 
how the JV can be recognized in the IP. 

For any feature with a size ~ Lg comparable with the 
healing length to become optically resolvable it must be 
enlarged (typically about 10 times) during the expan- 
sion. We consider the situation when axial expansion of 
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FIG. 3: Interference patterns of two expanded overlapping 
BEG clouds released from the waveguides as given by direct 
numerical integration of 13711 . Waveguides initially contained: 
(a) uniform BEGs, (b) two DS aligned at a; = 0, (c) JV solu- 
tion located at a; = 0. The image (d) is the relative intensity 
between (b) and (c); v = 0.01, t = 100, d = 2"=*/^ zo = 3. 
The unit of x and y is the coherence length U. 



the cloud of length L can be ignored. This imposes the 
limitation ^ L in the chosen units. Wc also ignore 
quasi-lD thermal fluctuations. In other words, the phase- 
correlation length Lip is taken larger than the system 
size L. In reality this is too strong of a requirement — 
the discussed feature can be seen when L^> Ls- 

Employing the ansatz (|35I36|I . the density after time t 
(such that R{t) ^ d, that is, t ^ cP) becomes 



n{x,y, z) 



( sin(Y)cos(^)^A"(a;i)+ (38) 



, cos(^)sin(:^)^'(xi))P 

where the overall factor is dropped, and tp'{x) = 
Vl + i^tanh(sa;) and ip"(x) = Qo/ cosh(sa::) are, respec- 
tively, the real and imaginary parts of the JV solution 
m given by s = Qq — ±\^1 — iv. The signs ± are 

due to two possible directions of the current circulation 
in the JV. It is instructive to consider three distinctive 
cases: i) two uniform condensates, which can be repro- 
duced from ea. pSfl by setting ip' — and -ip" — 1; ii) two 
identical DSs a,t x — Q, which can be obtained by setting 
-0" = 0; and iii) the JV. In the case i), taking the limit 
L — > cx) and performing explicit integration, one finds 



the well known parallel fringes n{x,y,z) ~ cos^(zzo/i). 
In the case ii), there is a zero at a; = 0. Its width can be 
estimated from ea. (l38|l as 5x y/i in the limit y/i ^ 1/s, 
that is, when the DS length Lg — 1/s has expanded sig- 
nificantly: Ls <^ t/Lg. The actual density profile can be 
obtained analytically for 1/s ^ < \/t as 



nDsix,y,z) 



4(l + i/)e" 



i3 



.:e2^os2(^). (39) 



It features the parallel fringes with the central zero as 
shown in Fig. 3b. 

The JV IP can easily be understood by analyz- 
ing the vicinity x = 0. Indeed, for 1/s <C 
l^^l ^ Vt, the tanh(sa;) function can be replaced by 
a step function and Qq / cosh(sx) effectively becomes 
S{x) J (IxQq/ cosh{sx) — t:Qq6{x) I s. Thus, the density 
profile due to the JV becomes 



njv{x,y,z) 



-(2vl + i'cos(— ^) X 



t3 



± 



■sin(— )) . 



t 



(40) 



The profile of the zeros of the density n = n{x, z,t) de- 
fines the feature specific for the JV. In the DS case, zeros 
belong to the set of mutually orthogonal lines a; = and 
z — 7r(rt + l/2)</zo, with n integer. In the JV case, rep- 
resented by ea. (|4()|l . the lines of zeros do not cross any 
more and obey the condition 



x = ±- 



/ 1 - 3t/ 



tan 



(f) 



(41) 



The inclined tangential feature seen on Fig. 3c is a con- 
sequence of the smooth relative phase change from to 
2tt in the JV. It is also worth noting that the JV circu- 
lations in different directions produce tangential slope of 
different sign in the IP. 

On Fig.© we have plotted column densities (inte- 
grated over y) at t — 100 (t — 50ms in usual units, which 
is a typical experimental expansion time after which the 
absorption image is taken). It should be noted that the 
above presented IPs correspond to the case when the sep- 
aration between the waveguides Zg is significantly larger 
than the transverse extension d. Obviously, in this situ- 
ation the tunneling v is essentially zero. If one tries to 
decrease zq /d, the visibility of the fringes worsens due to 
the exponential factor in H38|l so that for zo/d ~ 1 just 
one central fringe is seen. Thus, in order to achieve a 
good resolution, the clouds should be quickly separated 
from each other and, then, immediately released. Un- 
der these conditions, the JV solution formed at a closed 
proximity between the waveguides will have no time to be 
distorted by the inter-particle interactions after the tun- 
neling is cut off. Obviously, the duration of the waveg- 
uides' separation from some distance Zq ~ d (when the 
tunneling is finite) to, e.g., zq ~ lOd (when the tunneling 
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is, practically, zero) is limited from below by the inverse 
frequency « (P of the radial confinement. From above, 
it is limited by the axial response time w l'^. This 
requirement can safely be satisfied if {d/lcY *C l,that is, 
in the quasi-lD regime. 

The above discussion has been limited to destruc- 
tive imaging of the JV. Very recently, a non-destructive 
method has been employed 13], 14] to continuously sam- 
ple the relative phase of two spatially separated BECs. 
Tilted fringes in the interference pattern of the outcou- 
pled matter waves have been seen as an evidence of axial 
gradients of the relative phase 0] . This method can also 
be very usefuU for detecting the Josephson vortex and its 
conversion into the dark soliton and vice versa. The out- 
coupling pulse produces recoil atoms characterized by the 
wavefunction which, in the co-moving frame, is a replica 
(apart from the numerical factor) of the wavefunction of 
the confined atoms. Thus, the interference pattern pro- 
duced by the outcoupled clouds will be identical to the 
one discussed above. 



V. CONCLUSION 

The atomic Bose Josephson vortex can be created by 
the phase imprinting technique and detected due to its 



particular feature in the column density by absorption 
imaging performed after some ballistic expansion. The 
Josephson vortex can be controllably displaced by impos- 
ing tunneling current (created by disbalance of chemical 
potentials) between the waveguides. In quasi-lD, motion 
of an atomic Josephson vortex is strongly coupled to the 
current circulation through the phase-slip effect. This 
leads to a destruction of the circulation for the vortex 
speeds above a certain value determined by the Joseph- 
son coupling. In contrast to the standard approach to 
Josephson vortices in superconductors within the Sine- 
Gordon formalism, a description of the coupling between 
the center of mass motion and the circulation necessarily 
involves both density and phase variations. 
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